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ABSTRACT 

Two-particle annihilation reaction, A + A — > inert, for immobile reactants on the 
Bethe lattice is solved exactly for the initially random distribution. The process reaches 
an absorbing state in which no nearest-neighbor reactants are left. The approach of 
the concentration to the limiting value is exponential. The solution reproduces the 
known one-dimensional result which is further extended to the reaction A + B — > inert. 
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Recent studies of reaction-diffusion systems have emphasized fluctuation effects 
and breakdown of the standard chemical rate equations in low dimensions. For the 
simplest reactions of two-particle coagulation, A + A — >• A, and annihilation, A + A — > 
inert, on the one-dimensional lattice, several exact results have been reported [1-13]. In 
the diffusion-limited, instantaneous-reaction case, these processes show non-mean-field 
power-law decay of the A-particle density. 

Another solvable limit, in one dimension [14], is the case of no diffusion at all. 
Generally such models of immobile reactants have received less attention in the recent 
literature [14-16]. The reason is that unless longer-range reactions are allowed for [15- 
16], the time dependence involves exponential relaxation to an absorbing state rather 
than power-law behavior typical of the fast-diffusion reactions. Thus there are no 
universal fluctuation effects involved. 

On the other hand, immobile-reactant systems provide an example of freezing in an 
absorbing state with an explicitly nonuniversal, initial-condition-dependent behavior 
persistent at all times and, again, not consistent with the mean-field rate equations. 
It is therefore of interest to derive exact results whenever possible. Thus far, they 
were available only in one dimension [14]. In this work we report an exact solution 
for A + A — > inert on the Bethe lattice. We also derive exact results for the reaction 
A + B — > inert, limited to one dimension. 

Examination of the one-dimensional solution of A + A — > inert, [14], suggests 
close similarity to the models of random sequential adsorption [17-19]. Specifically, 
the annihilation reaction resulting in removal of two nearest-neighbor fixed reactants is 
equivalent to "deposition" of a "dimer" of two empty (reacted) sites. The two processes 
are dual to each other in that nearest-neighbor pairs of sites available for deposition 
correspond to unreacted pairs of reactants. Mathematically, the only difference is in 
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the initial conditions. In deposition, the lattice is usually assumed empty at time t = 
which would correspond to the full occupancy for reaction. 

This connection to random sequential adsorption models suggests that exact so- 
lutions can be sought for models of multiparticle annihilation reactions corresponding 
to n-mer deposition [19], and for models formulated on high- connectivity lattices such 
as the Bethe lattice, etc.; see [18,20-21]. In this work we consider the Bethe-lattice 
case; the methods of [20-21] are adapted for reactions. Besides notational differences, 
this essentially amounts to a more careful treatment of random initial conditions. The 
latter is achieved by a method different from the techniques used in deposition model 
studies [20-23]. Since the Bethe lattice of coordination number 2 is identical to the 
one-dimensional lattice, we obtain the d = 1 solution as well. However, further simpli- 
fication in d = 1 allows us to also solve exactly the reaction A + B — > inert. 

The Bethe lattice of coordination number z > 2 can be viewed as an interior 
part of the infinite Cayley tree: each site is connected by bonds to z nearest neighbor 
sites, and there are no closed loops formed by bonds. In fact the details of the lattice 
connectivity are not important for our considerations. However, we disregard any 
end-effects. Since the number of "boundary" sites in a finite-number-of-generations 
Cayley tree grows proportional to the total number of sites provided z > 2 (i.e., there 
is branching in each generation), the boundary effects can be profound when long- 
range spatial correlations are present such as at phase transitions, both static [24] and 
dynamical [25]. The model considered here does not have any "dangerous" spatial 
correlations; size effects can be safely ignored. 

Consider a fc-site connected cluster on the Bethe lattice. One interesting feature 
of a loopless lattice, shared with the d = 1 lattice for which z = 2, is that in such a 
cluster the k sites are connected by exactly k — 1 internal bonds. This statement is 
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well known and easily established by induction: each new site can only be connected to 
one existing cluster site, by one bond, because loops are not possible. Another useful 
conclusion is that the number of bonds shared by the cluster sites and the nearest 
neighbor sites immediately outside the cluster under consideration, is zk — 2(k — 1). 
Here zk is the total number of neighbors seen by all the k cluster sites, while the term 
2(k — 1) subtracts the number of neighbors internal to the cluster (twice the number 
of bonds). 

In the model of immobile reactants we assume that the lattice sites are initially 
occupied at random with probability p, and empty with probability 1 — p. The initial 
reactant density per site is c(0) = p, and we would like to calculate the time-dependence 
of the density c(t) at later times t > 0, given that each nearest-neighbor reactant pair 
annihilates with the rate R per unit time. In fact, it is convenient to absorb the rate 
in the dimensionless time variable 



Similar to the random sequential adsorption studies [17,22] and some recent results 
for particle-exchange dynamical models [26-27] , we consider the probability (r) that 
a connected fc-site cluster is fully occupied by reactants at time t. Configuration of 
the sites which are nearest-neighbor but exterior to the cluster can be arbitrary. Thus, 
initially, 



t = Rt . 



(1) 



P fc (0) = p" 



k 



(2) 



At times r > 0, the quantities -Pfc(r) remain the same for all cluster topologies 
on a loop less lattice, provided the initial conditions are topology independent. This 
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is because the topology dependence is not generated dynamically by the evolution 
equations for the Pfc(r). Indeed, their time variation is only determined by the number 
of possible reaction events within the cluster, which is equal the number of bonds, k — 1, 
and the number of possible reaction events in which the reacting pair has one site within 
the cluster and another outside the cluster. The latter again is not dependent on the 
cluster topology for loopless clusters. We argued earlier that the number of such pairs 
is zk — 2k + 2. Thus, the time-dependence can be obtained from the relations 



—^ = -(k- l)P k - (zk -2k + 2)P k+1 . (3) 

CIT 

Here the first term is self-explanatory. In the second term, the probability Pk+i is 
used because a larger cluster must be actually occupied in order for a reaction event 
involving a site outside the original fc-site cluster to proceed. 

Relations (3) apply for all k > 1. The most interesting quantity is Pi(r) = c(r). 
It is expected to decrease in time but remain finite as r — > oo. All other probabilities 
Pk>i are expected to vanish for large times. 

The solution of the recursion relations (3) can be obtained by various methods. 
Perhaps the simplest is to note that the Ansatz 



P k (r) = c(r)[a(r)] k - 1 , (4) 

where a(0) = p, eliminates the ^-dependence. Indeed, substitution in (3) shows that 
the solution of the form (4) is possible provided 



(5) 



dc 

— = —zac 
dr 



(6) 



The solution is 



a — 



1 + - 2)p (1 - e- T ) ' 



(7) 



c = 



p[l + (z-2)p(l-e- T )] 



z/(z-2) 



(8) 



For z = 2 the <i = 1 solution is obtained either as a limit of (8) or directly, 



2p(l-e- T ) 



(9) 



An interesting feature of expressions (8)- (9) is the explicit nonlinear dependence of the 
surviving reactant density for all times on the initial density p. The approach to the 
limiting density as r — > oo is exponential, ~ e~ T . 

We now turn to the strictly one-dimensional case. Additional restrictions on the 
cluster topology now allow solution of the two-species reaction A + B — > inert. Let us 
assume that initially the d = 1 lattice sites are occupied by reactant species A with 
probability a, and by species B with probability /3, where a + (3 < 1. The sites are 
empty with probability 1 — a — /3. Nearest-neighbor AB and BA pairs react with rate 
R; see (1). 

As before, we consider the probability that a £>site-long cluster is "fully reactive." 
Let Ak(r) denote the fraction of fc-site clusters which are fully filled, with reactants 
in the configuration ABAB..., i.e., the leftmost site is A and the sequence is fully 
alternating. Similarly, let -Bfc(r) denote the fraction of fc-site clusters of the type 
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BAB A . . ., with the leftmost site B and otherwise fully alternating order. The state 
of the neighbor sites outside the /e-cluster is not important in the definition of these 
probabilities which are conditioned only on the internal cluster configuration. 

Let us introduce the quantities 



w = — ^— , (10) 



p=VaP , (11) 
where the notation p in (11) will become clear later on. Then initially we have 

^-evenW = B fc-even(°) = , ( 12 ) 



^fc-odd(°)/« = %odd( )/^ = P ■ ( 13 ) 

Let us denote the concentrations, per site, of reactant species A and B by a(r) = 
A 1 (t) and b(r) = -Bi(r), respectively. We have 

o(0) = a , (14) 



6(0) = (3 , (15) 



a(r) — b(r) = a — j3 
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(16) 



The latter relation is obvious. Thus we only have to calculate the sum a + b. It proves 
useful to introduce the quantities 



P k (T) = [A k (r) + B k (r)}/2 , (17) 



so that the sum a + b is given by 2P\ . 

Now the probabilities A k and B k satisfy the relations 



-JL = -{k- l)A k - A k+1 - B k+1 , (18) 



-JL = -(h- l)B k - B k+1 - A k+1 . (19) 

Here the first term corresponds to internal pairs reacting, the second term corresponds 
to the rightmost site reacting with an occupied external site, while the third term in 
both relations corresponds to the leftmost site reacting "externally." 

For P k we get, by summing (18) and (19), 



— ^ = -(k-l)P k -2P k+1 , (20) 
clt 

which is in fact identical to the d = 1 random sequential adsorption recursion [17,19,22], 
as well as to the z = 2 variant of (3). However, the initial conditions are more compli- 
cated, 



^-even(O) = / , (21) 
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(22) 



The solution is obtained by methods similar to solving the Bethe-lattice recursions. 
We try forms similar to (4), 

P k-even = ^even(r) [S(r)]*" 1 , (23) 

^fc-odd = ^odd(^) PW]*" 1 • (24) 

These, when substituted in (20), yield coupled first-order differential equations which 
can be solved explicitly to give 



2P k = /-^-^-^ Up + u,)e- 2 "( 1 - e ~ T ) + (~l) k (P ~ oo)e 2 ^- e ^) 



(25) 



The reactant concentrations then follow as 



a(r) = l\ a -(3+(p + ^e' 2 ^ T ) - (p - tu)e 2p ( 1_e T ) 



(26) 



f3-a + (p + u})e- 2p ^- eT ) - (p - u)e 2p ( l - e ~ T ) 



(27) 



where the parameters are related via (10)-(11). 

As in the fast-diffusion case, which was not solved exactly even in d = 1 but only 
analyzed asymptotically [7,28], the functional form of the concentration is different 
depending on if the initial concentrations are equal or not, although the difference 



here is less spectacular. For a = (3 (= ui = p), the time-dependence of each of the 
species concentrations is identical to the d = 1 result (9). The time-dependence via 
the double-exponential expressions entering the general results (26) and (27) is also 
similar to (9) with the effective concentration given by the geometrical mean p\ see 
(11). However, for a ^ (3 the full time-dependence is more complicated than for the 
single-species reaction, involving both the double-exponential term and its inverse. 

In summary we presented exact solutions for two-particle annihilation reaction on 
the Bethe lattice. For the coordination number z > 2, results were reported for the 
single-species case. In d = 1, where z = 2, we reproduced the known exact single- 
species result and extended the solution to two-species reaction. 
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